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Abstract 

The Probabilistic Dynamic Synthesis method is a technique for obtaining the 
statistics of a desired response engineering quantity for a structure with non-deterministic 
parameters. The method uses measured data from modal testing of the structure as the 
input random variables, rather than more "primitive" quantities like geometry or material 
variation. This modal information is much more comprehensive and easily measured than 
the "primitive" information. The probabilistic analysis is carried out using either response 
surface reliability methods or Monte Carlo simulation. In previous work, the feasibility of 
the PDS method applied to a simple seven degree-of-freedom spring-mass system was 
verified. In this paper, extensive issues involved with applying the method to a realistic 
three-substructure system are examined, and free and forced response analyses are 
performed. The results from using the method are promising, especially when the lack of 
alternatives for obtaining quantitative output for probabilistic structures is considered. 
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Introduction 


Accounting for the statistical geometric and material variability of structures in 
analysis has been a topic of considerable research for the last 30 years. The determination 
of quantifiable measures of statistical probability of a desired response variable, such as 
natural frequency, maximum displacement, or stress, to replace experience-based "safety 
factors" has been a primary goal of this research. There are several problems associated 
with its satisfactory application to realistic structures, though. The first problem is 
accurate definition of the input random variables (rv’s). The random variability of the 
material characteristics is generally available, but the variation in the geometry of a complex 
structure, such as a hollowed-out turbine blade, is virtually impossible to define accurately. 
The second problem for analysis is the large size of the finite element models frequently 
used to simulate these structures. Substructuring reduction techniques are required to 
decrease the size of many structural systems just to perform a single, deterministic analysis. 
Finally, a significant problem in the analysis of probabilistic structures is the accurate 
generation of the Cumulative Distribution Function (CDF) necessary to obtain the 
probability of the desired response variable. Monte Carlo (MC) simulations can be 
performed to calculate probability distributions of the desired output, but several thousand 
runs are required for accurate results. 

This research applies a methodology previously developed by the authors, called 
"probabilistic dynamic synthesis" (PDS) to solve these problems 1 . The PDS method uses 
dynamic characteristics of substructures measured from modal test as the input random 
variables rather than "primitive" random variables such as material and geometric variability. 
These characteristics, which are the free-free eigenvalues, eigenvectors, and residual 
flexibility, are readily measured and for many substructures, a reasonable sample set of 
these measurements can be obtained. The variation in these dynamic attributes accurately 
accounts for the entire random variability of the substructure. This is in contrast to using 
“primitive” rv’s, which no matter how well measured, cannot capture the variation of every 
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geometric grid point location in addition to the variation in material properties. Using the 
residual flexibility method of component mode synthesis (CMS), these dynamic 
characteristics are used to generate sample models of the substructures that are coupled to 
form sample system models, which are considerably smaller than an unreduced finite 
element model. Finally, these sample models are used to obtain the CDF of the response 
variable by either applying Monte Carlo simulation or by generating datapoints for use in 
the response surface reliability method, which saves a substantial amount of computer time. 

The previous work by the authors verified the feasibility of the PDS method on a 
simple seven degree-of-freedom (dof) spring-mass system. In this paper, we examine the 
extensive issues involved with applying the method to a realistically-modeled three 
substructure system. Both free response and forced response analyses are performed. The 
goal of the free-response analysis is to quantify a +/- 3a band about the median natural 
frequency, compared to the somewhat arbitrary +/- 5% band about the deterministic 
solution frequently used in industry. The goal of the forced response analysis is to quantify 
a + 3a value for the maximum response of a dof on the structure, compared with the 
completely experience-based "factor of safety" presently applied to the maximum 
responding dof in a deterministic analysis. These results are presented by comparing a 
CDF obtained using a Monte Carlo simulation of the "baseline" model, in which "primitive" 
random variables are used, with a CDF using the PDS method. While there is considerable 
room for improvement, these results are promising, especially when existing options for 
design of realistic structures are considered. Finally, potential sources of error and avenues 
of further research are presented. 

Probabilistic Background 

The research described in this paper employs the response surface reliability method 
approach of probabilistic structural mechanics for determining the statistical structural 
response characteristics. A more extensive review of the reliability methods is presented 
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by Brown 2 . To briefly review this technique, consider a limit state function g(X) = Y(X) - 
y, where y is specific value. Cornell 3 , and Hasover & Lind 4 developed the First Order 
Reliability Method (FORM) by dividing the X space into two parts, g < O (Y < y) and g > 
O (Y > y), and approximating g as a first order Taylor series expanded about the mean of 
each rv. 

If X is transformed to independent standard normal rv's U, then this first order 
approximation allows the multi-dimensional probability density function (PDF) to be 
represented by a one dimensional Gaussian PDF (see Fig. 1). Therefore, P(Y< y) = 0(|3), 
where <&(.) is the Gaussian CDF function found in handbooks, and P is shortest distance 
from the g=0 curve to the origin, located at the Most Probable Point (MPP). This greatly 
simplifies the calculation of probability values for a given limit state. The reliability method 
was expanded by Rackwitz 5 to multi-dimensional problems for which the limit state curve 
g=0 is an explicit, nonlinear function of the rv's, which makes the determination of p much 
more difficult. The method makes use of Lagrange techniques in an iterative algorithm to 
find p. Rackwitz & Feissler 6 and Chen & Lind 7 continued the development of this method 
by examining how to develop an "equivalent normal distribution" for rv's with non-normal 
distributions. 

The application of FORM and its extensions for non-explicit limit state functions, as 
is the case for large structural finite element models, requires the use of numerical 
differentiation to obtain a first order approximation of the performance function Y(X). A 
CDF can be constructed from this first order approximation using the technique described 
above or directly from the first order approximations of the mean and standard deviation if 
the distribution is assumed to be Gaussian. A practical second order approach was 
developed in 1987 by Wu & Wirshing 8 , who improved the accuracy of probability levels 
obtained by FORM for non-explicit limit states by using a partial second order expansion 
called the Advanced FORM (AFORM). Initially, the FORM is performed and MPP's 
obtained. The limit state is then expanded about each MPP using only the first order and 
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the pure (no mixed variable terms) second order terms of the Taylor series expansion. The 
first and second order terms are then "linearized" using a change of variable, and the MPP 
search algorithm used in FORM is applied. A CDF can also be generated directly from this 
partial second order approximation if a lognormal distribution is assumed. Wu also 
developed the Advanced Mean Value (AMV) method, a procedure for updating the FORM 
or AFORM solution by using the original, exact solution for g. All of these reliability 
methods have been implemented into a new probabilistic finite-element computer code 
developed by Southwest Research Institute, NESSUS 9 . 

System Description 

In order to examine the applicability of the PDS method for a realistic design 
problem, a structural system had to be chosen that met several criteria. First, the structure 
had to be modeled using standard methodology. This was satisfied by using the widely 
used commercial finite element code, NASTRAN. Next, since tying CMS together with 
probabilistic methods is an important goal of this research, the size of the model had to be 
large enough to be able to realize a substantial reduction in computer time due to dynamic 
reduction. This requirement is satisfied by using a structural system composed of a "disk," 
which is made up of 630 quad4 plate elements and constrained at the center, and two 
"blades," which are each composed of 24 quad4 plate elements (see Fig. 2). To pave the 
way for future mistuned bladed-disk analyses, the disk was assumed to be deterministic and 
the blades non-deterministic. For ease of analysis, the standard assumption that in-plane 
translation and rotations are small has been applied. Finally, the structure had to possess a 
"generic" level of randomness not easily defined by a single variation in a material or 
geometric property. A single rv would not be able to capture variations in mode shapes that 
are independent of variations of natural frequency, for instance. To achieve this goal, each 
blade was separated into three sections, with two of the sections having a thickness set to be 
an independent normally distributed primitive random variable. In addition, the density of 
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each blade across all three sections was defined as a normally distributed independent rv, 
thus giving each blade three independent rv's (see Table 1). 

Monte Carlo Baseline Analysis 

A Monte Carlo analysis of the original, unsubstructured system using the assumed 
primitive random variable set was chosen to be a baseline for comparison with PDS. A 
1000 sample case was executed, taking over 20 hours of wall-clock time. In addition, MC 
analysis was required to simulate the modal testing phase of the PDS methodology to 
obtain the statistics of the dynamic rv's of the probabilistic blade substructures. The 
NESSUS computer code has a new interface with NASTRAN that manages the simulation, 
both by creating the random vector set for the given input random variables, and by 
automatically submitting the jobs. The input random variables can be any defined 
NASTRAN geometric or material property. 

Probabilistic Dynamic Synthesis 

The PDS methodology makes use of the residual flexibility method of CMS. This 
method has been developed by MacNeil 10 , Craig & Chang 11 , and Martinez et al. 12 . The 
essential idea in CMS is that substructure modes are truncated since their higher modes will 
not have a major effect on the combined system modes. The residual flexibility method 
incorporates the effects of the higher modes by determining their flexibility. A side benefit 
is that all the elements of the system stiffness matrix can be obtained from test and that the 
mass matrix can be closely approximated by a unity matrix in the non-boundary partition. 
Since all the information can be obtained from test, all the probabilistic information can be 
incorporated into the system matrices for response analysis. 

The first step of the probabilistic dynamic synthesis (PDS) method is to divide the 
model of a structure into substructures m= a,b,..,p, and the degrees of freedom into internal 
and boundary locations. Each substructure is represented by n samples, each of which is 



modally tested in a free interface condition. For this research, this testing step is simulated 
by a MC analysis of the structure using the defined primitive rv's. For substructure m, 
sample i, the test will yield eigenvalues and eigenvectors [Op 1 . The eigenvectors 

have to be consistent from one sample to the next, requiring an interfacing FORTRAN code 
to assign a consistent sign to the mode shapes for each sample and obtaining a consistent 
set of rigid-body modes (parallel to the coordinate system axes) for the two free-free modes, 
which is achieved by using the SUPORT card with the modified Householder method in 
NASTRAN. In addition, the boundary partition of the residual flexibility matrix [G bb ] mj is 
obtained from the measured boundary drive point frequency response functions of the 
boundary coordinates 13 . Here, the residual flexibility is calculated analytically by 
subtracting the flexibility for the retained modes from the total system flexibility (inverse of 
the stiffness matrix) for the constrained substructure. For the free-free substructures, the 
system is unconstrained so the stiffness matrix is singular. The "inertia-relief method,” 
initially developed by Craig 14 , is therefore used to obtain the residual flexibility. 

For the free-response solution using the PDS method, only the kept (non- truncated) 
eigenvalues, the boundary coordinates of the kept eigenvectors, and the boundary partition 
of the residual flexibility matrix are needed. These values can be combined into a single 
vector {x} mi , defined as 

{x} m ‘ = [x, ...x k {<t> b }T ...G bb ]"" T (1) 

where k is the number of kept modes. If the entire sample of substructure m is tested, {x} m 
can therefore be defined as a vector composed of elements that are each a random variable 
with measured mean and standard deviation. This vector is now transformed to {x'} m , a 
vector of standard normally distributed rv's, using the measured mean and standard 
deviation of each rv. An important assumption is made that the original rv's in {x} are 
normally distributed. This distribution is required for decorrelation of the random 
variables (discussed below), and is somewhat accurate, as shown in distribution-matching 
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tests published in the previous work by the authors 1 . If the measured distribution were 
found to match another standard type better, then an could be obtained. Since the match 
with a Gaussian distribution is not perfect, though, there is some error introduced. A 
technique to address this error will be discussed later in this paper. Another method of 
increasing the accuracy would be to obtain an “equivalent normal” distribution using using 
the Chen-Lind three- parameter method 7 In addition, there will be some degree of 
correlation between each of the random variables, which can be calculated from the 
measured data. This information is placed in a correlation matrix [C] m relating each element 
with every other element. For the FORM analysis, a set of independent random variables 
{ u } m is required. This can be accomplished by making an orthogonal transformation of 
{xT with the eigenvectors of the correlation matrix to uncouple the {x’} coordinates, 
thereby creating {u} m . This can be expressed for substructures m =a,b,..,p as 

{x’} m = [0] c m {u} m (2) 

It becomes evident at this time that the size of the dynamic rv set is intractable for a 
realistic problem. For this case, which has only two probabilistic substructures and where 
20 modes are retained (out of a possible 108) per substructure, the number of dynamic rv's 
is 802. Several assumptions are therefore made to drastically reduce the number of 
dynamic rv's. The first is to assume that the limit state function is insensitive to the variation 
in the rotational dof s in the modes. This reduces the size of eigenvector rv's from 240 to 80 
per substructure. The second is to assume the limit state is insensitive to not only the 
rotational dof s, but also to the off-diagonal terms in the boundary residual flexibility matrix, 
which reduces the size of that contribution from 144 to 4 per substructure. These 
assumptions do not remove these variables from the formulation of the substructure 
stiffness matrices; instead, it allows the use of the mean value of those variables (or median, 
as will be discussed later). Admire, et. al 15 examined the effect of completely removing the 
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off-diagonal G elements, and calculated a natural frequency error of less than 5% between 
the exact value and the value obtained using the residual flexibility formulation. 

In addition, a substantial reduction in the number of independent rv’s is achieved by 
eliminating "artificial" rv's. These were first found by calculating the correlation matrix of 
the independent rv's [C] u . Since the elements of [C] u are uncorrelated, it was expected to be 
an identity matrix. The result was somewhat different; for the rows/columns of [C] u that 
corresponded to eigenvalues on the order of one, the diagonal term was equal to 1.0 and the 
other elements in the associated row/column were very small. However, many of the 
eigenvalues turned out to be very small and the analogous rows and columns in [C] u were 
not diagonal; the diagonal values were 1.0 but the other terms in the row/column were all on 
the order of one (between zero and one) instead of being close to zero. This can be 
explained by realizing that the actual number of independent rv's in the system is equal to 
the original number of independent primitive rv's, not the number of dynamic rv's. 

Therefore, the rv's associated with the small X, c values are actually insignificant and the 
magnitudes of these values are mainly due to numerical roundoff. The high off-diagonal 
values in the new correlation matrix in these rows/columns therefore indicates a large 
amount of correlation between these numerical artifact "rv's". The above result can be used 
to further decrease the size of the rv set {u}. A sum of the eigenvalues of the original 
correlation matrix is calculated, and any eigenvalue less than 3% of the sum is deemed 
insignificant and the associated rv in {u} neglected. The value of 3% was reached by 
starting at 5% and decreasing the cutoff value (thereby including more rv's) until the final 
result stabilized. This reduced the eigenvalue set for each blade from 101 to 1 1. It is unclear 
why this value did not actually reduce to three, which is the number of primitive rv's per 
blade, but nevertheless, this reduction greatly facilitated the analysis. 

To generate a first order Taylor series representation of the limit state, each 
independent random variable is varied individually by some percentage of its standard 
deviation a, which is equal to the square root of the corresponding eigenvalue of the 




correlation matrix, while the other rv's are kept constant at their mean values. For higher 
order methods, each independent rv is varied by two different amounts to allow a quadratic 
curve fit. Matrix [u] is created by varying each of its elements by +/- 0.3 of a standard 
deviation, which for that rv is the square root of the correlation matrix eigenvalue. The 
number of rows in [u] is equal to the number of probabilistic substructures times the 
reduced number of independent rv’s, and the number of columns is equal to this number 
times two plus one additional column for the mean value set, as shown below. The 
superscripts on the rv u refer to the number of the probabilistic substructure, and the 
subscript refers to the number of the independent rv. 


[u] = 


-3c , 

“i 

0 

0 


3c , 

u l 

0 

0 


0 0 0 0 
0 0 0 
0 -.3c . .3c . 0 

uj uj 


(3) 


Each column in [u] is then transformed to the set of correlated standard normal rv’s {x} 1 
using the transpose of [<1>] C , and then into the original rv’s {x} using the standard normal 
transformation. The new vectors {?*.}, [<J>], and [G bb ] are pulled out from {x} and placed in 
substructure mass and stiffness matrices according to the residual flexibility formulation: 


+ d> kT G _1 <t> -d> kT G _1 

T VJ bb M/ bk ' J bb 


K“ = 


M = 


sym 


’bb 


'I O' 
0 0 


(4) 


where A k is a diagonal matrix of the kept eigenvalues {A.}\ For the purposes of this study, 
the third natural frequency of the combined system was selected as the response value of 
interest, so the system matrix is then compiled and an eigensolution performed to obtain this 
value for each column of [u]. The columns of [u] along with the associated response value 
are then input into the NESSUS code where the AFORM partial quadratic method can be 
applied. By simply changing the first variation in the above approximation to a value very 
close to zero times the standard deviation, the same program can also be easily altered to 
create a FORM linear approximation. 
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After initial examination of the results showed a poor comparison with the MC 
baseline, an important subtlety was found in the application of the transformation to 
standard normal variables. For the MC baseline solution, the results are obtained by 
creating a normal distribution about the mean values of the primitive rv’s; this distribution is 
symmetric, so the median is equal to the mean. Since the finite element solution is slightly 
nonlinear as a function of the input rv's, though, the mean of the solution will not equal the 
median, but the median solution will result from using the mean (equal to median) primitive 
rv's. 

As PDS is implemented for these analytical simulations, though, there is an 
intermediate step that introduces error. Symmetric Gaussian MC distributions are created 
about the means of the dynamic rv's, but these rv's are actually not purely Gaussian, since a 
slightly nonlinear eigensolution has been performed to obtain them, and the means will not 
equal the medians. This results in a skew of the entire CDF curve since the method assumes 
that the dynamic rv's do in fact follow a symmetric distribution. The error can be 
substantially corrected by simply ensuring that the median primitive values are carried 
through the procedure to generate the median results. This is accomplished by using the 
medians x of the calculated dynamic rv's as the "mean" in the transformation to standard 
normal coordinates: 

(x)={x’}o+ x (5) 

This change drastically improved the results in both the mean and standard deviation of the 
results since it moves the expansion of the limit state about a more accurate location in the 
u-space. 

A parallel application of the PDS method is to apply the dynamic rv's using the MC 
approach rather than the reliability approach. The MC should provide more accurate 
answers, while the reliability method will be more computationally efficient. This technique 
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is performed by using the statistical data obtained using the modal testing simulation to 
create normally distributed samples for each column in [u] rather than creating "load cases" 
for the Taylor series. This will result in a [u] matrix with the number of columns equal to 
the number of Monte Carlo samples. 


Results 

The PDS FORM and PDS MC results for the CDF of the system third natural 
frequency are compared with the MC baseline in Figure 3 and Table 2. This frequency was 
chosen because the associated two nodal diameter mode shape is similar to the type of 
excitation shape found in turbomachinery which the forced response analysis would 
simulate. The higher order reliability methods were also applied to this problem, but they 
did not show any improvement over the linear case; this may be due to the very small 
nonlinearity of the limit state. From a design problem approach, the errors as seen on the 
table for the .01 and .99 CDF levels are less than 2%, so the values could be used with 
confidence. From a theoretical standpoint, the curves do not line up extremely well; this is 
indicated by the error in the standard deviation of about 20%. Potential sources of error 
include the following: 1 ) transformation error of the statistics of the original primitive rv's 
due to inaccuracy of the MC simulation of the dynamic rv's (possibly insufficient number 
of samples), 2) transformation error in the statistics and recorrelation of the rv's due to 
their assumed normality, 3) error due to truncation of the number of independent rv's, 4) a 
small truncation error due to using the RF formulation of CMS (examined for the 
deterministic case), and 5) programming error. 

It is interesting to compare computational times at this point. All the analyses were 
performed on a CRAY T-90 supercomputer. The MC baseline analysis took 2527 CPU 
seconds and 20.3 hours of wallclock time to run. PDS MC took 389.5 seconds of CPU 
time and 1 1 12.3 seconds (18.5 minutes) wallclock time, a drastic reduction enabled by the 
continuous generation, storage, and analysis of modally reduced system matrices rather 
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than the complete model generation and analysis necessary for each of the MC baseline 
samples. The FORM method took only 20.1 seconds of CPU and 64.2 seconds of 
wallclock to create the datapoints for input to NESSUS, and about 30 seconds wallclock to 
run NESSUS to perform the FORM solution. This further reduction is due to the 
substantial decrease in the number of solutions generated. As anticipated, both the use of 
CMS and probabilistic methods drastically decreased the amount of time necessary to run 
the analysis, which, especially for a more detailed model, could be a pre-requisite for use in 
design. 


Forced Response 

The use of PDS was now expanded to forced response. In particular, a frequency 
response solution of the system was derived since this type of excitation would prove most 
applicable to the bladed-disk problem. For forced response, the dofs of a substructure are 
partitioned into three sections, x 0 , which are internal dofs with no external load applied, x,, 
internal dofs with an external load applied, and x b , boundary dofs which may or may not 
have external load applied. The residual flexibility transformation matrix from the original 
set of coordinate to the generalized set is therefore: 


X. 

1 


<D.-G ,G 1 
i res — ob res- 

0 

0 


bb^b 


0 

<t>. -G , G _1 
l res - lb res - bb b 

0 


i .G ^ , , 

res — ob res-bb 


G 


-1 


=[T] { q } (6) 


res - ib res - bb 

I 


The transformed mass and stiffness matrices will be identical to those obtained for the free 
response shown in Equations (4). To determine the right-hand side of the forced response 
equation of motion, the transformation matrix and the original load vector are both carried 
throughout the analysis until the generalized force is needed, at which time the matrix 
multiplication is carried out. Therefore, in addition to the boundary partitions of the modal 
matrix and the residual flexibility matrix required for the free response solution, the 
partitions associated with the dofs with load applied are required for the forced response 
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analysis. To more easily perform this analysis, the transformed coupled system is 
transformed to an uncoupled system using the set of truncated mass-normalized system 
eigenvectors [O], resulting from the free response solution of the system, as shown below: 

{q },=[<*>] >{q] 2 • (7) 

The resulting uncoupled system is in the form 

[I]{q} 2 + [C]{q} 2 + [A], {q} 2 = K [T] T {f } (8) 

where [A]j is the diagonal matrix of the eigenvalues of the transformed system. The 
standard (viscous) assumption of constant modal damping of 0.5% is used for [C], The 
input force is a system of four harmonic loads, each of amplitude 1 lb. (4.448 N) and 
applied in-phase with a peak displacement of the two nodal diameter mode shape to simulate 
bladed-disk excitation. There are standard techniques available for obtaining the frequency 
response solution for this set of uncoupled single-dof systems. Because of the series of 
linear transformations performed on the solution, though, a complex quantity for the 
solution is desired rather than one in terms of magnitude and phase, as is commonly derived 
in texts on the subject and which cannot be easily transformed. This complex vector {q} 2 
can now be transformed back to the original coordinates {x} using equations (6) and (7): 

{x}=[T][<D] 1 {q} 2 (9) 

The absolute values of this complex vector are then calculated to obtain the physical 
response values. 

As seen above, applying PDS for the forced response solution requires more 
information from the modal data and the residual flexibility matrices. This increases the 
number of dynamic rv's necessary to solve the problem, which increases the complexity of 
the problem in several aspects. To minimize this number, it is important to decide a-priori 
which internal dof s will either have external load applied or require a displacement recovery. 
Although the entire modal matrix and residual flexibility matrix are calculated, only those 
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partitions of the modal matrix and the residual flexibility matrix, along with the boundary 
dofs, are stored. In addition, only the correlations with these additional dofs are generated. 

To address the design problem defined earlier, the response variable for the forced 
response analysis was chosen to be the maximum response for all dofs in the structure. 
Initially, the excitation was applied at only a single frequency, but this was changed to a 
wider excitation bandwidth since the maximum response occurs at the damped natural 
frequency of the mode shape being excited, and this mode shape can occur over a range of 
frequencies for a probabilistic structure. It is assumed that the actual excitation mechanism 
could also vary in frequency by this amount, which is generally the case in engine 
turbomachinery, for instance. Using the maximum over some frequency range also helps to 
reduce the extreme variability in response that can be observed at any particular excitation 
frequency 16 , which would introduce substantial nonlinearities in the response limit state 
surface. A range of +/- 10 hz about the deterministic natural frequency of 458 hz was used 
in this case. This value is probably inadequate, since the free response case shows that the 
actual range is between 438 hz and 472 hz, so this may be a source of error in the final 
results. 

The Monte Carlo baseline frequency response solution required extensive alteration 
of the MSCPOST subroutine within the NESSUS code. The PDS procedure using both 
the MC and reliability approaches for the frequency response problem is similar to the free 
response problem. The main additions are the subroutines necessary for calculation of the 
different [T] matrices for each substructure and the solution algorithm for the complex 
frequency response. Because the value of the largest responding dof is no longer at a 
frequency known a-priori (since the third natural frequency and third natural mode will vary 
for each statistical sample), a sorting algorithm to find the maximum value was used that 
scanned the response values for all the selected dofs for all the frequencies in the bandwidth 
chosen. The dofs selected were those at the blade tips, which were assumed to contain the 
maximum responding dof of the entire structure for the chosen excitation. This maximum 


16 



displacement value was then used as the NESSUS response variable, and input along with 
the corresponding load set. PDS MC and FORM cases were run as well as a normal 
distribution using the linear approximate mean and standard deviation. The PDS quadratic 
AFORM method using the medians was also attempted, but the AFORM algorithm in 
NESSUS/FPI never converged to a solution, which is a problem that sometimes occurs with 
the technique 17 . 

The CDF results for the PDS methods are plotted along with the 1000 sample MC 
baseline simulation in Figure 4 and Table 3. The MC baseline CDF is extremely 
unsymmetric, so this clearly causes a problem for the reliability methods. This asymmetry 
may be due to the fact that the maximum response is obtained by scanning the response of 
different nodes, and this may skew the distribution. Further examination of this phenomena 
was performed on a more detailed bladed-disk, which will be reported in the future, and 
additional work is also needed. It is also noted that there is fair agreement between the MC 
baseline result and that obtained from MC using dynamic rv's. The errors for the design 
point at .99 are still around 10% for both PDS methods, though, which is within reason 
especially when one considers that at present, there are no alternative methods for obtaining 
this maximum value. As with the free-response analysis, the significant error in the 
standard deviation and the lack of coincidence of the curves indicated visually signifies that 
there is still error present in the methods. An additional source of error in the reliability 
methods for this analysis is the extent to which the maximum responding dof does not have 
a smooth dependence on the rv's. Possible discontinuities in this dependence would cause 
errors in the response surface generation created using numerical differentiation. An 
examination of the run times yields similar results as the free response case, with the MC 
baseline taking 27.3 hours of wall-clock time, the PDS MC taking 18.5 minutes wall-clock, 
and the PDS FORM taking about 90 seconds. Because of the acceptable level of error for 
design, the increased accuracy in system variability obtained by using dynamic rather than 
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primitive rv’s, and the tremendous reduction in computer time, the technique was applied to 
the solution of a mistuned bladed-disk, which will be reported in a future paper. 

Conclusions 

A probabilistic dynamic synthesis method has been applied to the analysis of a 
realistic three-substructure system. The new technique uses modal data from a sample set 
of substructures to generate a set of dynamic random variables which fully describe the 
probabilistic variation in the structures. The residual flexibility method of component mode 
synthesis is then used to generate probabilistic mass and stiffness matrices which can be 
used to obtain any desired response variable. Probabilistic analysis is performed on these 
stochastic systems by applying both the traditional Monte Carlo technique and new 
reliability methods. Solutions for both free and forced response of the realistic system were 
obtained, and the results compare favorably with a baseline Monte Carlo analysis. The new 
method results in a dramatic decrease in computer run-time as well as increased accuracy in 
representation of the probabilistic variation of the parameters of the structural system. 
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Figure Captions Page 

Fig. 1 Joint probability density surface 
Fig. 2 3-Substructure system 

Fig. 3 3-Substructure mode 3 natural frequency comparing Monte Carlo with PDS 
Fig. 4 Percent errors for forced response solution 
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Table 1 3-Substructure System Information 


Deterministic 
Number Elements 
Number Nodes 
Number dofs 
Number unconstrained 
dofs 

Boundary Conditions 

Material 

Geometry 


Fixed at center and two 
adjacent nodes 
Steel, E= 206. 8 MPa, 

v=.29, p=7.85 g/cm 3 ) 
Diameter = 6.375 cm, 
Thickness t = .254 cm 


Blade B 

N on-deterministic 
Number Elements 
Number Nodes 
Number dofs 
Boundary Conditions 

Material 

Geometry 


N on-deterministic 
Number Nodes 

24 

Baseline Svstem 

Number dofs 

36 

Number Elements 

Boundary Conditions 

108 

Number Nodes 

Material 

Free-Free 

Number dofs 

Geometry 

Steel, E= 206.8 Mpa 

Number uncon. Dofs 


v=. 29 

Boundary Conditions 


p ~ N(7.85 g/cm 3 ,.785 
g/cm 3 ) - independent rv 
Length =5.08 cm 

Others properties 


Thickness: 

same as above; all rv’s 


Section 2: t = .254 cm 

independent 


Section 3: 

t - N(. 254, .0254 cm) - 
independent rv 
Section 4: 

t ~ N(.254,.0254 cm) - 
independent rv 



24 

36 

108 

Free-Free 

Steel, same properties as 
Blade A 

Length = 5.08 cm 
Section 5: t = .508 cm 
Section 6: 

t ~ N(.508,.0508 cm)- 
independent rv 
Section 7: 

t - N(.508,.0508 cm) - 
independent rv 


678 

695 

2085 

2076 

Fixed at center and two 
adjacent nodes off center, 
node 2 and 12 





